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The recoil effects of spontaneous photon emissions during optical pumping of a trapped three-level atom are 
exactly calculated. Without resort to the Lamb-Dicke approximation, and considering arbitrary detuning and 
saturation of the pump laser, the density of recoil shifts in phase space is derived. It is shown that this density 
is not of Gaussian shape, and that it becomes isotropic in phase space only for a branching ratio corresponding 
to fluorescence scattering but unfavorable for optical pumping. The dependence of its anisotropy on the laser 
saturation is discussed in the resonant case, and the mapping of moments of the atom's center-of-mass motion 
due to the pumping is presented. Moreover, it is shown how optimum parameters for protecting the center-of- 
mass quantum state from pump-induced disturbance depend on the specific property to be protected. 
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I. INTRODUCTION 



Optical pumping (OP) QJJ, has proved to be a very power- 
ful method for preparing specific electronic states of atoms. 
By application of a near resonant light field the atoms are 
excited and subsequently decay into the desired electronic 
state via spontaneous emission. In bulk media such as a 
solid, or a gas of atoms or molecules, optical pumping can 
be understood solely in terms of electronic or optical proper- 
ties of the medium. Although in dense gases modifications 
arise due to atom-atom interactions mediated by the radia- 
tion field, the mechanical effects of optical pumping may be 
safely neglected in such situations, since momentum is trans- 
ferred to the collective sample. However, any light interac- 
tion with free or weakly bound atoms not only is accompa- 
nied by a change of the internal electronic quantum state of 
the atoms, but unavoidably also affects their center-of-mass 
motion lHHD]. Examples of this type of interaction are found 
with the momentum diffusion in fluorescence scattering of 
free atoms J5L S_ 3, la Bal , with the laser cooling of trapped 
atoms fiol fill Il2 Jl3l \\M or ions lfl5l flrjll . The mechanical 
effects are particularly apparent when we consider the inter- 
action of a laser field with an individual particle, e.g., with a 
single trap ped ion lfl7ll . with single atoms confined in a mi- 
cro cavity Il8l [l9(|. as well as in atom optics l20ll and laser- 
control of atoms [210. When OP is applied to such a system, 
one has explicitly to take into account the mechanical effects 
of the laser-atom interaction. Sometimes, as with OP of low- 
lying vibronic states of weakly bound trapped ions, the radia- 
tive recoil generates a desired effect: cooling of the sample. 
However, more often than not, the mechanical effect modifies 
the momentum or energy distribution of the sample and the 
concomitant radiative spectrum in an undesirable way. In par- 
ticular, the incoherent nature of spontaneous emission makes 
it add decoherence I22I1 to the quantum state of the center-of- 
mass motion of the atoms. In view of the possible applica- 
tion of such systems to quantum-state engineering l23l l24ll . 
quantum logic, and computing 125l Ba, |27ll . the introduction 
of any such decoherence is obviously problematic. Thus, the 
identification and avoidance of decoherent effects is of utmost 



importance. 

From the experimental viewpoint, the optimum parameters 
must be found in order to obtain smallest decoherence. Given 
a certain atomic species confined in a trapping potential, sev- 
eral parameters are predetermined such as the electronic tran- 
sition frequencies, the relaxation rates, and the trap frequency. 
Apart from modifying the trap, the adjustable parameters are 
the laser intensity and detuning. One may consider conve- 
nient to use fast optical pumping, faster than the vibrational 
period, since the negligible center-of-mass motion improves 
the efficiency of the pump process. On the other hand, the re- 
quired laser intensity may saturate and broaden the employed 
atomic transition, which spoils OP's resolution of vibronic 
states. Thus, the interplay of saturation and pumping time 
requires close inspection. 

The recoil upon atoms interacting with light has been 
treated in various contexts, such as for example in the light 
interaction of free |9] and trap ped sing le atoms IT28ll or in the 
collective atomic recoil laser 1291 l30fl . The applied models 
consider an atomic two-level system continuously irradiated 
by near-resonant monochromatic light. Such a treatment of 
free atoms leads to a Gaussian momentum distribution in ac- 
cordance with the central limit theorem applied to the statis- 
tically independent individual photon recoils [8]. For trapped 
atoms this situation is different since here the waiting-time 
distribution of subsequent spontaneous photon emission cor- 
relates the individual photon recoils in phase space, and the 
validity of the central-limit theorem is doubtful. Additionally, 
OP requires a third atomic level, a lower target state. The 
process is intrinsically terminated by the spontaneous decay 
into this state. Thus, contributions with a small number of 
radiative events may dominate the interaction. The recoil ef- 
fects of these contributions make the resulting distribution of 
the atom in vibrational phase space substantially deviate from 
being Gaussian - the more so, the less is the number of inter- 
action events with the pump light. 

In this paper, we address a single atom trapped in a har- 
monic potential. We model the motion of the atom's center-of- 
mass by a ID quantum-mechanical oscillation along a certain 
direction. This model can be viewed as a prototype system 
that provides us with the essential information to be applied 
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to one of the real systems mentioned above. Using a quantum 
trajectory approach we obtain analytic results for the quan- 
tum statistical properties of the atomic center-of-mass motion. 
Moreover, we obtain insight into the decoherence properties 
by employing a phase-spatial picture, that lends to illustrative 
interpretation. 

The method of quantum trajectories l3ll [32j [33j [34j [35j[36j 
l37ll is perfectly matched to studying OP with proper inclu- 
sion of recoil effects, and avoiding the Lamb-Dicke approxi- 
mation: Firstly, the trail of spontaneously emitted photons of 
a single experimental realization is finite, since a stationary 
state is reached when the system is no longer affected by the 
laser drive. This fact simplifies the structure of single quantum 
trajectories. Secondly, when using waiting-time distributions, 
easy connection is made between the center-of-mass vibration 
and the statistics of photon emissions. Finally, an elegant for- 
mulation of recoil effects is found that characterizes the pro- 
cess of optical pumping, from the viewpoint of the vibrational 
motion, by just a single distribution function. 

In Sec. |II]we introduce the model system under considera- 
tion and develop its evolution as a quantum trajectory in phase 
space. Section Hn] outlines the vibronic state of the pumped 
atom in terms of the density of recoil shifts in phase space. 
This density is characterized in Sec. [IV] by its statistical mo- 
ments that are expressed by quantities derived from photon- 
counting statistics. Section [V] is then devoted to a discussion 
of the results in limiting cases, such as fluorescence scattering 
and maximally anisotropic scattering in phase space, as well 
as to optimum laser parameters in order to minimize detri- 
mental recoil effects of OP. Section [VI] places the results in 
the context of possible observations. Finally, Sec. I VIII pro- 
vides a summary and conclusions. 



II. QUANTUM TRAJECTORIES OF OPTICAL PUMPING 
A. Quantum master equation 

Our generic model for treating the motional effects of OP 
consists of a A-type three-level atom whose center-of-mass 
coordinate is bound in a harmonic trap potential. In this way, 
photonic -recoil effects can be described for various systems, 
such as single ions in rf traps or neutral atoms in optical-dipole 
or magneto-optical traps. The level scheme and setup of opti- 
cal pumping is shown in Fig. Q] 

In both equivalent schemes, a laser field is assumed to drive 
the electronic dipole transition |1) <-> |3), whereupon photons 
are spontaneously emitted connected with electronic dipole 
relaxation from the auxiliary level |3) to both |1) and |2). The 
latter process occurs only once, since in state |2) the atom de- 
couples from the driving laser and thus has reached is final 
electronic state. Both level schemes (a) and (b) include pump- 
ing of population from the electronic state |1) to state |2), 
with the corresponding direct transition being dipole forbid- 
den. The optical-pumping rate is determined by the electronic 
relaxation rates 71 and 72 of the two dipole transitions, cf. 
Fig. Q] and by the laser Rabi frequency k = —2dE/h. Here, 
d is the electronic dipole moment of the transition |1) <-> |3), 
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Figure 1 : Draft excitation schemes for the inversion of the two-level 
system |1) — > |2) by optical pumping via the third auxiliary level |3). 
Scheme (b) is identical to (a) with the ground state being replaced by 
a metastable state. 



and E is the electric-field amplitude of the laser. 

The subject to be addressed is the atom's center-of-mass 
vibration at frequency v along one principal axis of the trap- 
ping potential. Furthermore, only the motional recoil effects 
of spontaneous photon emissions are of interest here. The ad- 
ditional recoil effects due to the laser interaction have been 
discussed elsewhere |38ll and can be separated by choosing 
the propagation direction of the laser beam perpendicular to 
the direction of motion under consideration. This motion may 
be thought to extend along the axis of a linear trap, while the 
orthogonal confinement of the particle be "strong", such that 
the trap absorbs the momentum transmitted by the laser. 

The dynamics of the vibronic density operator g(t) is given 
by the quantum master equation ll38ll 



do 
di 



±.[H Q + V(t),g} 
in 



(1) 



J2 ^ 



where {^4, B} = AB + BA denotes the anti-commutator, and 
Ho describes the free time evolution of the atom in the trap 
potential JH], 



H a = hvPb 



0=1,2 



2ftlO a (Ja,z 



(2) 



with uj a (a = 1, 2) being the electronic transition frequencies. 
Moreover, b and ¥ are the annihilation and creation operators 
for the vibration along the chosen principal axis of the trap, 
and the electronic spin-type operators are given by 



a a - = \a}(3\, <5a,+ = 
*a,, = I (|3><3| - |a><o|) . 



(3) 
(4) 



Finally, the external drive of the laser with frequency u> is de- 
scribed by the Hamiltonian 



V(t) 



1 



h.c. 



(5) 



with k = 2d3i • E/fi, being the Rabi frequency, CI31 = (3|d|l) 
being the transition dipole moment, and E being the electric - 
field amplitude of the laser. 
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In Eq. ([T]i the super-operators 7Z a are responsible for the 
momentum recoil due to spontaneous emission of a photon 
via the electronic decay channels |3) — ► \a) (a = 1, 2). They 
act on an arbitrary operator A as 



TZ a A= I ds^ a {s)e ls ^ x Ae 



(a=l,2) (6) 



with the respective photon-momentum recoil hk a = hui a /c 
(a = 1,2) being projected on the considered axis according to 
the normalized dipole radiation characteristics, 



/■2-7T 

■ / deb {1 - [n a • n(ft)] 2 } , (7) 
Jo 



where s = cos 8 is the projection of the photon-emission di- 
rection on the considered motional axis. The unit-normalized 
vectors n a and n(fi) point into the directions of the transi- 
tion dipole moment (3|d|a) and of the spontaneous emission 
il = (8, <fi), respectively. Assuming, for example, transition 
dipole moments orthogonal to the direction of motion under 
consideration, Eq. Q reduces to 



*,(-) = J (l + «") 



(8) 



The strength of the recoil effect of individual spontaneous 
photon emissions can be characterized by Lamb-Dicke pa- 
rameters: The position operator x is expressed in terms of 
the annihilation and creation operators as k a x = rj a (b + b>) 
(a = 1,2). The Lamb-Dicke parameters r\ a of the corre- 
sponding transitions |3) — ► \a) are defined as 



i] a = k a Ax (a =1,2). 



(9) 



where Axq = y/h/(2mv) is the rms spread of the position 
of the atom of mass m in the ground state |0) of the trapping 
potential. 

To proceed, the master equation is transformed into a frame, 
whose center-of-mass part oscillates with the trap frequency v 
and whose electronic part rotates with the laser frequency to. 
The transformed master equation is written in terms of three 
super operators, 



dg 
di 



M + Ji{t) + J 2 {t) 



Q, 



(10) 



where the evolution in the absence of spontaneous photon 
emissions is described by 

AfQ=^_(H eS Q-QHi s ). (11) 
The occurring non-Hermitean effective Hamiltonian reads 



H eS = -(«ffi,+ + K*a-i,-)+fiA|l)(l| 

-— 2^ laVa, + Va,-- (12) 



with A = oj — u)\ being the detuning of the laser. 

The jump operators J a {t) describe spontaneous photon 
emissions via the transitions |3) — > \a) (a = 1,2) and are 
given by 

Ja(t)g = J a S a ,-n a (t)g, (i = l,2) (13) 



where the spin super-operator reads 



+ • 



(14) 



and the now time-dependent recoil operators read 



n a {t)g= / dsn a (s)V(i Va se lut )g. 



The action of the rhs super operator is defined as 



(15) 



T>(a)g = D{a)gD*<(a), 



with the coherent displacement operator, 



D(a) = exp(aa^ — a*a). 



(16) 



(17) 



Thus, in Eq. ( fT5l ) the length of the shift in motional phase 
space depends on both the Lamb-Dicke parameter rj a and the 
direction of photon emission s, whereas the phase of the shift 
is determined by ir/2 + vt, i.e. it depends on the time of the 
spontaneous emission event. 



B. Formal solution of the propagator 

The formal solution of the master equation ( fTQb in the ro- 
tating frame reads 



g{t) = Moi 2 (t,t )g{t ), 



(18) 



where g(to) is the initial density operator, and the propagator 
is given by the time-ordered expression 



M i2(t,t )=T exp jy dr [Af + Ji(t) + J 2 (t) 

(19) 

Eq. dT9b can be rewritten as a sum over all possible trajectories 
where the spontaneous emissions |3) — > |2), i.e. the operator 
J 2 (t), occurs at the random times ti, . . . , t n 15^1 : 
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■Moi2(*,*o) = ^2 dtn 



n=0 ° z ° 



dt 1 M -L(t,t n )J 2 (t n )Moi{t n ,t n -i) . ..J 2 (ti)Moi(ti,t ). 



(20) 



to 



r 



The operators M.oi(t, t') determine the time evolution be- 
tween the spontaneous emissions |3) — > |2) and contain both 
the laser drive and the spontaneous emissions |3) — > |1): 



M 01 (t,t') =Tcxpl / dr 7V r +Ji(r 



(21) 



After a spontaneous emission |3) — > \2) at a time 
tk, the atom is in its final (pumped) electronic state 
1 2) where it decouples from the laser drive, and thus 
M i(t k+1 ,t k )J 2 {t k ) = J 2 {t k ). Since J 2 {t k+1 )J 2 (t k ) = 
a second spontaneous emission is impossible — it would re- 
quire the electronic state to be |3) . Accordingly the propagator 
d2"0t reduces to only two parts, 



M i2{t,t ) = .Moi(Mo) 



dt'j 2 tf)M 01 (t',to). (22) 



The first term on the rhs in Eq. d22l corresponds to the quan- 
tum trajectory when no spontaneous transition from state |3) 
to 1 2) occurred. The second term represents the trajectory with 
one final spontaneous transition |3) — * |2). Only one such 
spontaneous transition can take place. 



Defining now the projection operators into the electronic 
states la), 



(23) 



and using the facts that (i) the projection into the pumped state 
1 2) is not affected by the dynamics on the transition |3) «-> 1 1), 
that is V 2 M i(t,t') = V 2 , and that (ii) V 2 J 2 {t') = J 2 (t'), 
the part of the propagator taking the atom to the pumped elec- 
tronic state 1 2) is obtained fromEq. (1221 as 



V 2 Moi 2 (t,t )=V 2 + dt'J 2 (t')Moi(t',to). (24) 

Jto 

The first term represents the initial population in level |2), and 
the second one the effect of the pumping. 

Analogous to Eq. d20l i the propagator M 01 (t,t') can be 
formally rewritten as a sum of sequences of spontaneous pro- 
cesses 1 3) — ► Denoting by t n = t' the time of the final 
spontaneous emission that leads into the electronic state |2), 
Eq. (1241 thus can be rewritten as 



00 rt nt 2 

V 2 M 012 {t,t ) =V 2 + J2 dtn--- dt 1 J 2 (t n )M (t n ,t n . 1 )J 1 (t n . 1 )...A(h)M (ti,t ). (25) 
n =\Jta Jto 



Here the operator A4o(t, t') describes the evolution exclusive 
of any spontaneous emission, 

M (t, t') = cxp \£[{t - t') \ . (26) 

A non-vanishing action of Ji{t{) requires population in 
level |3). Since A^o(ii,^o) does only couple levels |1) and 



1 3), the second term in Eq. d25l l only acts on populations in 
levels |1) and |3). Assuming zero initial population of g(to) in 
the rapidly decaying state |3), the action is restricted to popu- 
lation in level |1), and we may thus multiply the second term 
in Eq. (l25T l on its rhs with the projector V\, 



V 2 M 012 



d*lJ2(*n)A<o(*n,tn-l)Jl(t n -l) • • • M) (*1, to)V X . 



r 



(27) 



The second term describes sequences of n consecutive n excitation processes from level |1) to level |3), as shown 
spontaneous transitions. The laser interaction is restricted to by the sequences of operators M Q alternating with jump op- 
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erators J x oc S\ ,_ = |1)(3| . . . |3)(1|, see Eqs ©, ©, and 
(TBI . Each of these excitation processes results therefore in 
the waiting-time distribution Q S3 EH E] 



w(t-t') = ( 7 i+72)(3| M (t,t')\l)(l\ |3) 



(71+72) (3\U eS (t-t')\l) 



(28) 



obtained as the probability (density) of the atom being in level 
1 3) at time i, on the condition that it was in |1) at time t', and 
that no photons have been emitted in the interval \t' , t]. Here 
we have defined U e s(t) = exp(—iH e gt/h). The waiting- 
time distribution is calculated for resonant pumping in App. 
ID1 Using definition (f28b the formal solution < f27b can therefore 
be written as 

00 ft rt 2 

V 2 M ai2 (t,t ) = V 2 + J2 / dt n--- / Mi 

n -±Jt J to 

x[x 2 n 2 (t n )w(t n - f n _i)] 



x |[Ai7ti(t n _i)iy(f n _i - t n - 2 )] ■ ■ ■ 
...[Ai^i^Mti-to)]}^! (29) 



where 



A a =7o/(7i+72), (a = 1,2) 



(30) 



are the branching ratios, and the operator that induces the tran- 
sition into the final level |2) is defined as 



V21Q = {&2,-0~l,+) Q (0-1,-0-2,+) 



(31) 



Defining the joint probability density for n spontaneous 
transitions at times t\, . . . , t n as 

w n {t n , . . . ,h) = A 2 A"~ 1 u; (tn - t n -i) . . . (32) 
...wfa-t^wfa-to), (n>0) 

the result (f29l is rewritten as 



°° rt 

V 2 M i 2 {t,t ) = V 2 + Y^ / dt n ... I dt 1 w n (t n ,...,t 1 )K2(t n )Tl 1 (t n - 1 )...K 1 (t 1 )V2i 

n=l 



(33) 



m. THE VIBRATIONAL EFFECTS OF THE OPTICAL 
TRANSITIONS 



where the total recoil shift is 



In Eq. d33l . a sequence of recoil operators, cf. Eq. ( fT5l ), 
appears, which can be combined as 



TZ 2 (t n ) . . .TZl(ti) = j dSnfl 2 (Sn) ■ ■ ■ J dsifii(si) 

x P(n ?2S „e Il/t ")P^ lS „_ 1 e Il/t "- 1 ) 
...P(i?/isie wtl ) 

ds n [i 2 (s n ) ...I dsin\(s\) 



i({t n , s n }, {ti, si}) = ir] 2 s n e ll/tn + ^2 iViSme 1 



(35) 



The n-fold integration represents the averaging over the direc- 
tions of recoil in the n consecutive spontaneous emissions at 
times ti, . . . , t n . In order to consider any number of sponta- 
neous emissions at arbitrary moments of time within the in- 
terval [to, t], we must average over the joint photon-emission 
probability at the consecutive emission times. In this way we 
xT> [a n ({t n , s n }, . . . , {ti, s\})] , (34) obtain the operator of the mean recoil shift in phase space: 

I 



°° rt rt 2 

V(t,t ) = J2 *»•••/ dt 1 w n (t n ,...,t 1 )Ti 2 (t n )TZ 1 (t n ^ 1 )...n 1 (t 1 ) 

n= \Jta J t 

= ^2 dtn dSn 

A*2 ( s n )■■■ dti ds 1 fi 1 (si)w n (t n , . . . ,ti) 

n= \ J to J J to J 

XV [a n ({t n , S n }, . . . , {tl, Si})] . 

I 



(36) 



One may define the density of recoil shifts in phase-space 
that contains all the recoil effects of the spontaneous emis- 



sions at random times into random directions, 

00 ft /> r-t 2 f 

p(a;t,t ) = / dt n J ds n ... / dt\ / ds\ 
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x M2(s n )Mi(«n-i) ■ ■ ■ fi>i(si)w n (t n , ...,h) 
x 8[a-a n ({t n ,s n }, . . . , {ti, si})] . (37) 



The operator of the mean recoil shift in phase space ( [361 ) is 
expressed as an integral over T> weighted by this distribution 
p(a; t, to) m phase-space: 



T>(t,t ) = / d 2 ap(a;t,t )T>(a). 



(38) 



Since p(a;t,to) is an integral over positive probabilities, it 
may be viewed as a blurring distribution that smoothes the 
initial vibrational phase-space distribution. It contains all the 
characteristics of the random process that acts on the center- 
of-mass degree of freedom of the electronic population that is 
transferred from state |1) to state |2). 

In the propagator of Eq. (133) . leading to the pumped elec- 
tronic state 1 2), the second term on the rhs contains the aver- 
aged displacement as given in Eq. ( [36) . and thus Eq. ( |38l ) may 
be employed to yield 



V 2 M Q u(t,t ) =V 2 + / d 2 ap(a;t,to)V(a)r: 



'21- 



(39) 



Applying Eq. d39b on the initial vibronic density operator 
g(to) and tracing over the electronic degree of freedom, we 
obtain the reduced motional density operator for the popula- 
tion in the pumped electronic state, 

(2|e(t)|2) = (2|e(to)|2>+ / d 2 ap(a;t,t )t>(a)(l\g(t Q )\l). 



(40) 

A vibrational density operator (a|g|a), that corresponds to the 
population in electronic level |o), is related to the s-ordered 
phase-space distribution 



^i s) («) = ^Tr 



6 {s \a-b)(a\g\a) 



(a = 1,2), 



(41) 

where the s-ordered delta operator is formally defined via the 
displacement operator (TT7T > as l43ll44ll 

6 {s) (a-b) = J d 2 f3exp(a(3* - f3a* +s\(3\ 2 /2)D(l3), 

(42) 

and the probability to find the atom in the electronic state \a) 
is defined by 



P a = Tr((a\g\a)). 



(43) 



We want to finally come up with a mapping of the initial 
vibrational phase-space distribution on the distribution gener- 
ated by the optical pumping process. For this purpose we ap- 
ply Eq. (|4TT > on the mapping of vibrational density operators, 
Eq. ( l40i i. making use of the definition ( TToT ) and obtain 

P ( 2 \a 7 t) = P 2 (t )P [ 2 s) (a,t ) + [ d 2 a'p(a';t,t Q ) (44) 



xTr \^ s) {a - b)D(a'){l\g(t )\l)D\a') 



Employing the cyclic property of the trace, one rewrites the 
trace as 



Tr 



D^a')S^(a-b)D(a')(l\g(t )\l) 



= Tr 



a-(b + a') (l|e(t )|l) 



(45) 



where the displacement operation results in a shift of the vi- 
brational operator b — > b + a' . This result, with substitution 
of integration variable a' — > a — a' in Eq. d44l) . yields 



P 2 {3) (a,t) = P 2 (t )P^ s> (a,t )+ I d 2 a'p(a- a';t,t ) 



xTr 



(l|ff(*o)|l) 



(46) 



With the use of definition ( f4Tb this phase-space distribution 
writes 



P 2 (s) (a,t) = iM*o)iT(Mo) 



(*)/ 



(47) 



+P 1 (t ) / d 2 a'p(a-a';t,t )p( s \a' 1 t ). 



Thus the initial phase-space distribution of the population in 
electronic level 1 1) is mapped onto a contribution to the phase- 
space distribution of the pumped level |2) via a convolution 
with the density of recoil shifts in phase-space (|37| |. 



IV. THE PHASE SPACE DENSITY OF RECOIL SHIFTS 



A. Moments of the density of recoil shifts 

The density p of recoil shifts in phase space multiplied 
by an infinitesimal phase-space volume, p(a; t, to)d 2 a, is the 
probability of displacing, by a, the initial phase-space distri- 
bution within the time interval [to,t]. The observable features 
of the collective recoil shift will be characterized by this den- 
sity's statistical moments of all orders (k, I), 

{a* k a l ) p{a . tM ) = J d 2 aa* k a l p(a;t,t ). (48) 

To evaluate these moments we may employ the characteris- 
tic function p(a, t, to), which is the Fourier transform of the 
density of recoil shifts, 



p(a;t,t )= d z pexp(a(3* - 0a*)p(/3;t,t o ). (49) 



With the help of this characteristic function the statistical mo- 
ments, Eq. d48l i. are obtained as j45ll 



(a* k a l ) p(a . t , to) = (-1)' [d*& a .p(a, t, t )] q=0 . (50) 

Inserting Eq. d37b into d49l ) and performing the integral over 
(3 yields 
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p(a;t,t ) = jr f dt n ( ds n ... ( * dh ( d Sl (51) * Y[ cxp {^ai m s m e - c.c.) . 



exp [~mr] 2 s n e 

n-l 
m—1 



x fJ-2(s n )^i(s n -i) ■ ■ ■ fJ-(s{)w n (t n , ...,ti) 
x exp [aa* n {{t ni s„}, . . . , {ti, si}) - c.c] , 

where a„({i„, s„}, . . . ,{ti, si}) is defined in Eq. ([35). With 
this definition, the exponential factor writes 



exp [aa* {{t n , s„}, . . . , {ii, si}) 



c.c. 



cxp 



Each exponential factor on the rhs can be rewritten as 

exp {—air\se~ lvt — c.c.) 

= cxp [-2ir 1 s\a\$t(e wt - laisa )] 

= exp [— 2i?ys|a| cos(«4 — arga)] , 

so that Eq. ( IBTt reads 



p(a;t,t ) = y^ dt 



n=l 
X 



to 



dhw n (t n ,...,h) / ds n p 2 (s n )e- 2 ^ s ^ cos ^"-^^ 



ds„_i//i(s„_i)e^ 2i,)lS "- l|Q|cos(l/t "- 1 - al ' SQ) ... / d Sl /ii(si)e- 2 "' lSl|Q|cos(l/tl - arsQ) 



(52) 



(53) 



(54) 



Each integral over a dipole characteristics is a directional 
spectrum of recoils; it is formally expressed by the Fourier 
transform 



MO 



dsii a (s) exp(-is£). 



(55) 



In Eq. d54l l this directional spectrum appears n times with 
the arguments 277 a |a| cos(vt m — arga), where m — 1, . . . , n, 
and a = 1 , 2 depending on the decay channel. Each of these 
spectra defines a characteristic function associated to the in- 
stantaneous probability density for a single recoil shift at time 
t, 



p (a,t) = p [2r) a \a\cos(vt — arga)] , (a =1,2) 



With these replacements Eq. d54l i takes on the form 



°° pt r-t-2 

p(a;t,t ) = I dt n ... d^w^t 

1 J tn J tn 



(56) 



,fi) 



n=l J z 



xp 2 (a,t„)£ 1 (o!,i„_i) . . .^(a-jti). (57) 

The moments of the density of recoil shifts are evaluated 
according to Eq. d50l . The characteristic function d57l i de- 
pends on a via the factors p (a, t). The derivatives of these 
factors (App. [A} are 

'd k a d l a .p a (a,t)] = (-i^^+'e^C-*) ((cos9) k + l ) a 

(58) 

with a = 1,2. The moments of the projection of the sponta- 
neously emitted photon's wave vector on the dipole axis are 

({cos0) k ) a = f dss k fi a (s). (59) 



Since the dipole radiation characteristics shows even spatial 
symmetry, p a (—s) = p a (s), cf. Eq. dA4K all odd moments 
vanish, ((cos 9) 2k+1 ) a = (k = 0, 1, 2, . . .). In consequence 
all odd-order derivatives of p ([3, t) vanish, see Eq. (1581) . and 
via Eq. d5"7l i the odd-order moments of p vanish likewise: 



*k„l\ 



}p(a;t,t ) 



= if k + l = 2n+ 1. 



(60) 



This result could have been anticipated, since the incoher- 
ent scattering process connected with optical pumping cannot 
generate a coherent vibrational amplitude. 

The non-vanishing even moments of a complete optical- 
pumping process (t — > 00) are 



(a* k a l 



= ( _1 )']C/ dtn... dt lWn (t n , . . . ,ii) 

x \d k d l a , p 2 {a,t n )p i {a,t n - 1 )...p i {a,t 1 ) } a _ Q - ( 61 ) 

For the lowest-order non-trivial even moment (a*a) p two 
derivatives have to be performed on any same factor p (a,t) 
which results in a sum of n terms, where the first one reads 

- ^ / dt n ... dtiw n (t n , . . . ,ii) (62) 
n —iJt J t 

x I \d a d a *p 2 (a,t n ) £ 1 (a,i„_i)...£ 1 (a,ii)j . 

Given thatp a (0,t) = /ujO) = 1 [Eqs <f35} and d56]l] and 

[d a d a *p Ja,t n )] a=a = —t;|((cos6*) 2 )2, the above term is 



7 ?2 2 ((C0S#) 2 ) 2 £ / dt r , 
n=l Jt ° 



dtiW n (t 7 



to 
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= r 1 l((cos9) 2 ) 2 Y J Pn 

n=l 

= V 2 2 ((COB0) 2 ) 2 . 



where to = was set to cancel the trivial phase factor e 2wt °, 
and the spectral waiting-time distribution is defined as 



(63) 



w(ll>) 



Here the statistics of the number of spontaneously emit- 
ted photons during the pump process is introduced as time- 
ordered integral of the joint probability density of photon 
emissions, 



dtw{t)e i 



(70) 



with w(t) defined in Eq. 

The last complex-valued factor in Eq. 



Pn 



dt n ■ 



to 



dtiw n (t„,...,h), (n> 1), (64) 



\ 2 w(2is) 
1 - \iw(2v) 



Ae l 



(71) 



to 



and Pq = 0. The moments of this photon statistics are evalu- 
ated in App. [B] 

The sum of the other n — 1 terms is 



- where A is a real-valued semi-positive number - depends 
on the branching ratio, and, via the waiting-time distribution, 
on the trap frequency and the laser parameters. This func- 
tion characterizes the saturation of the recoil distribution's 
anisotropy. With the definition 



E 



dt, 



dtiW n (t r , 



Th — 1 



w(2v) = wexp(i0 2ii ), 



(72) 



P =i 



and since w is semi-positive, the modulus A and the phase 
<p 1 (a,t n - 1 )...\d a d a . Pl (a,t p )]...p 1 (a,t 1 )} . (65) of the anisotropy (ED can be written as 

L — 1 J — 1 ) a— 

\ 2 W 

As the differentiated factor is — rj 2 {(cos9) 2 )i in all n — 1 
terms, they sum up to yield 



A = 



^((cosfl) 2 )!^ / dt n . 

n=l ^ to 

oo 

! ((cos0) 2 >i^P„(n-l) 



dt\w n (t n , . . . ,t\)(n - 1) 



tan < 



Employing Eqs 



(1 - Aiw) 2 + 4Aiwsin 2 (<^/2) 
sin 6 W 



is 4>w_ — XiW 

and d7Tb . we rewrite the moment 



(73) 
(74) 



-n p Ae i4,A . 



= ?7 f((cos0) 2 )i((n ph >-l): 



(66) -pjjg phase-dependent quadrature is defined as 



where the mean number of emitted photons is given in Eq. 
(IB41 i. Thus, the lowest-order non-trivial moment of the den- 
sity of recoil shifts is 



q(<j>) = -j= (ae'* + a*e~ 



14>\ 



as 
(75) 

(76) 



(a a) 



ii{(cos6) 2 ) 2 + r ? 2 ((n ph ) - l)((cos0) 5 



(67) 



Upon insertion of the mean number of spontaneously emitted 
photons, see ( 1B4I ). this moment becomes 

n p = (a*a) p = r 1 2 ((cos6) 2 ) 2 + ^7? 2 ((cos 0) 2 )i, (68) 



where the mean recoil-induced excitation n p = {a*a) p is in- 
troduced. This moment is the width of the density distribution 
of recoil shifts to which the system has been subjected in the 
optical pumping process. 

In the same order one obtains the moment (a 2 ) p which 
gives information on the rotational symmetry of the density 
of recoil shifts in phase space. Analogously to the above cal- 
culation one obtains this moment - see App. [C]- as 



The above considerations suggest that the average quadrature 
vanishes, (q((f>)) p = 0. Its variance, however, is in general 
non-vanishing and phase dependent. It is expressed in terms 
of the moments [Aq(<j>) = q(<p) — (q((j>)} p ] as 



([A^)f) p = («H + »((aV1, 

which, upon insertion of Eqs. d68l ) and d75l ). writes 

Sq 2 ^ = ([Aq^)] 2 ) p = n p [1 - A cos (20 + </> A )] 



(77) 



(78) 



where the abbreviation 8q p ^ = \J ([Aq((f))] 2 ) p has been in- 
troduced. 

At the values of the phase 0_ = —0^/2 and <f> + = <fi- + 
7r/2, the minimum and maximum variances, respectively, are 
attained, 



(« 2 ) P 



77 2 ((cos0) 2 ) 2 + ^ 2 ((cos< 

\ 2 w(2v) 
1 - Aiw(2i/)' 



n p (l ± A), 



and the product of the unequal rms uncertainties is 



(69) 



Sq P ,<p+ ■ Sqp.j,^ = n v \/\ - A 2 . 



(79) 



(80) 
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In general the density of recoil shifts is not rotationally sym- 
metric, as shows the ratio of the difference of maximum and 
minimum quadrature fluctuations over the mean fluctuation, 
i.e. the anisotropy in the quadratures: 



Sq 2 



Sq 



P,4>+ 



Sq 2 



A. 



(81) 



p,4>- 



There are, however, two special cases where the density is 
approximately isotropic, i.e. where A « 0. This situation 
is obtained either with negligible spectral waiting-time dis- 
tribution at twice the trap frequency, w = \w{2v)\ « 0, or 
with a very large value of the branching ratio, Ai « 1 and 
A2 ~ 0, see Eq. ( 1731 . The first case represents a situation 
where the spectral waiting-time distribution lacks a frequency 
component at twice the trap frequency. This condition means 
that the waiting times of the probabilistic sequence of sponta- 
neous emissions of photons are not synchronized to the half 
period of the trap oscillation. The latter case simply corre- 
sponds to a situation where the decay rate to state |1) much 
exceeds that to state |2), so that the optical pumping process 
requires a very long sequence of spontaneous photon emis- 
sions, each one randomizing the density of recoil shifts that 
eventually becomes isotropic. This case corresponds to fluo- 
rescence scattering, as discussed for a free atom in Refs |0,Ht]. 
In general, however, a finite sequence of photon emissions 
during the optical pumping generates a "squashed" density of 
recoil shifts with a predefined orientation in phase space that 
is determined by the phase of the spectral waiting-time distri- 
bution. 

Finally, for the characterization of the density of recoil 
shifts, the variance 



r 2 / *2 2\ -2 

dn p = (a a ) p - n p , 



(82) 



that corresponds to the mean n p , is evaluated. The first mo- 
ment on the rhs is obtained as [see Eq. ( IC12H 

(a* 2 a\ = ^((cos^ + ^acosfl) 4 )^ (83) 
+ 2n p (n p - ■ql{{cos6) 2 ) 2 ) (1 + Acas<f) A ) , 
so that the variance (l82l becomes 



Sn 2 



^((cos^ + ^tcostf) 4 ) 
A 2 



2n p (n p - T]l((cos9) 2 ) 2 ) (1 + A cos <f> A ) . (84) 



B. Mapping the moments of the atomic phase-space 
distribution 



So far, the density of recoil shifts has been characterized by 
its statistical moments. To characterize the vibrational quan- 
tum state of the ion after completion of optical pumping, the 
initial vibrational phase-space distribution of the ion has to 
be convolved with the density of recoil shifts, as given in Eq. 
(07}. Fourier transforming Eq. (l47l i. one obtains the mapping 
of the characteristic function, that corresponds to the vibra- 
tional phase-space distribution, as 

P 2 s \a,t) = P2(t )P 2 s \a,t )+P 1 (t )p(a;t,t )p!f ) (a,t ), 

(85) 

where the Fourier transform of the density of recoil shifts is 
given by Eq. S57\ . 

From the s-ordered phase-space distributions, Pa S \a,t), 
the s-ordered quantum-statistical moments are represented by 
integrals JHH, 



({# k (0(t)}) a = ^Tr [{bn i } s {a\g(t)\a 



d 2 aa* k a l p( s \ ai t), 



(86) 



where { b' b } s denotes the s-ordered product of operators and 
a = 1,2 specifies the electronic level. As a result of the 
Fourier transform, these moments are obtained equivalently 
as derivatives of the corresponding characteristics function, 



Q=0 



(87) 



Inserting the mapping of characteristic functions (1851 1 into 
( l87| i and applying the Leibniz formula yields the mapping of 
initial quantum-statistical moments of the ion's vibration on 
the final moments (to = 0, and t — > 00): 



({St fe (oo)6V)} ) =P 2 (0)({6t*( )6'(0)} ,) 9 + ^i(0)EE ( k r)(tr)( a<k ~ n)al ~ m ) ({^(°)^( )} } ■ 

S 2 " 2 n=0m=0 ^ ' ^ ' P " 1 

, m 

The moments of the density of recoil shifts, (a <x) p , act here Setting k = and I = 1 in Eq. d88l . the mapping of the 
as weight factors for contributions to the sum of initial vibra- coherent vibrational amplitude is obtained as 
tional moments on the rhs. 

(6(oo)>a = (6(0)), (89) 
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where the complete quantum-statistical average on the rhs is 
defined as 



= J2 Pa(---)a=Tv[g. 



(90) 



Thus the coherent vibration is unaffected by the pump process, 
as is the expectation value of the phase-dependent quadrature 
operator 



which reveals the mapping 



(90(°o)) 2 = (^(0)}. 



(91) 



(92) 



Although optical pumping is a highly incoherent process, the 
initial coherent amplitude of the ion's oscillation in the trap- 
ping potential is perfectly preserved. 

The mean vibrational excitation, however, is altered. Using 
Eq. < T88l > with values k = I = 1 and specifying normally 
ordered operator products, i.e. setting s = — 1, the mapping 
of the mean vibrational excitation is obtained as 



<ft(oo)} 3 = <»(0)}+Pi(0)np, 



(93) 



where n = b'b is the number of vibrational quanta and the 
mean number of vibrational quanta added by the pump pro- 
cess, n p , is given by Eq. ((68). Thus, as seen from Eq. d68l l, the 
addition of mean vibrational excitation by the optical pump 
process does entirely depend on the Lamb-Dicke parameters 
and on the mean number of spontaneously emitted photons, 
that is determined by the branching ratio, (n p h) = A1/A2 + 1. 
It does not depend on the speed of optical pumping that is 
determined by the parameters of the pump laser. 



Let us now turn to the rms spreads of these proper- 
ties. Whereas the expectation value of the phase-dependent 
quadrature is preserved during the optical pumping, see Eq. 
d92l ), the corresponding rms spread is not. The square of this 
rms spread, i.e. the variance, is constructed as 



(b 2 



(94) 



where only the last two terms are phase dependent. From 
Eq. d88l the mapping of the moment (b ) is obtained with 
the choice k = 0, 1 = 2 as 



(95) 



6 2 (oo)) =(b 2 (0))+P 1 (0)(a 2 ) p 



Thus, with the mappings d92l and d93l ), and after completion 
of pumping, the variance of the quadrature in the pumped 
electronic state becomes 



([A^(oc)] 2 ) 2 = <[Ag (O)] 2 > + Pi(O)^ 2 , 



(96) 



where the additional noise Sq 2 + is given in Eq. ( f78l . In accor- 
dance with the moments of the density of recoil shifts, mini- 
mum noise is added to the rms spread of the quadrature at the 
phase (/>_ = —<j) a I '2, and maximum noise is added at phase 
(j)+ = (j)- + tt/2. 

Finally most relevant is how broad grows the final vibra- 
tional number distribution. The variance of the vibrational 
quantum number is obtained from normally ordered moments 
as 



([An] 2 ) = (St 2 S 2 )-(n)((n)-l). 



(97) 



From Eq. ( I88I 1, using k = I = 2 and s = — 1, the mapping of 
the required normally-ordered moment is obtained as 



;St 2 (oo)6 2 M ) 2 = (^( )6 2 (0)) + Px(0) {25R [(« 2 ) p (& t2 (0))J + 4 (a*a) p (n^ + (a* 2 a 2 ) p ) 



(98) 



Thus using Eqs d93l and d98l > the mapping of the variance of also in its electronic ground state [-Pi(O) = 1, Pzfo) = 0], 



the vibrational quantum number, Eq. d97| l, results as 

([An(oo)] 2 ) 2 = ([An(0)] 2 ) (99) 
+ P 1 (0){5nl + 2n p [m 1 +P 2 (0)m 2 }} . 

The numbers mi and to 2 depend on the initial state, 



mi 



(«(0)>i + 5 



-A3? 



(6 2 (0)) 



ie 



-iipA 



m 2 = -f + (n(0)h-(mh, 



, (100) 
(101) 



has mi = 1/2 and m 2 = n p /2, the latter not entering the 
final variance. The mapping of the variance of the vibrational 
excitation then simplifies to 



([An(oo)] 2 ) 2 = n p + 8n 2 , 



(102) 



and Sri 2 , is given by Eq. 

For example, an atom being initially laser-cooled to its trap i.e., the variance of the final vibrational distribution equals the 
ground state [(n(0)}i i2 = (n(0)) = 0, (b (0)) = 0], and mean plus variance of the recoil displacements. 
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V. RESULTS 



A. Fluorescence scattering 

As noted before the case of fluorescence scattering is cov- 
ered by the model in the limit A2 — > 0, when the trail of 
spontaneously emitted photons becomes infinite so that the 
directions of recoils in vibrational phase space are efficiently 
averaged to result in an isotropic recoil density. However, the 
resulting recoil density would correspond to the limit t — ► 00, 
when the diffusion of fluorescence scattering would have in- 
finitely broadened the density distribution. Therefore, a direct 
comparison is impossible, and we restrict ourselves to a dis- 
cussion of the asymptotic behavior for small A2. 

In the limit A2 — > the anisotropy parameter becomes 

A-{ J' ^=1^^=°. (103 ) 

Apart from the exceptional case w(2u) = 1 IBUl . where the 
waiting times are perfectly synchronized to half the trap pe- 
riod so that a highly directional scattering occurs (A — ► 1), 
the recoil density is isotropic with respect to the variance of 
its phase-dependent quadrature (A — > 0). In more detail, 
for small A2 the anisotropy's asymptotic behavior is given by 
A x aA2 and </>a x const, where the coefficient a is 

w 

a= , = (104) 

yj(l - w) 2 +4wsin 2 (0 2i ,/2) 

The asymptotic form of the rms spread of the quadrature 
becomes then phase independent, cf. Eq. ( |78] l, and results as 

$q P ,4> X \f^p- (105) 

Like the mean vibrational excitation d68l i. the rms spread of 
the quadrature is now independent of the laser parameters. As 
the mean vibrational excitation asymptotically behaves as 




(106) 



A 2 



the rms spread of the quadrature ( 1105t asymptotically takes 
on the form, 

W((cos0) 2 )i/A 2 . (107) 

Correspondingly, the rms spread is asymptotically Sn p x n p 
and also turns independent of the laser parameters. Thus, 
we expect an isotropic recoil density in phase space that just 
weakly depends on the laser parameters through higher-order 
moments. 

Numerical evaluation of Eq. ( |37] i for t — ► 00 (to = 0) 
and A2 = 10~ 5 in resonant excitation (A = 0) reveals 
the isotropic recoil density in phase space as shown in Fig. 
12 which is independent of the preset laser saturation S = 
(\ K \/l) 2 > m agreement with previous considerations. Given 
the infinite trail of spontaneously emitted photons, each pho- 
ton contributing to a random recoil, one might presume the 




position 



Figure 2: Recoil density p(a) in dependence on position and mo- 
mentum scaled as x/(rjiAxo) and p/(r)iApo), respectively, with 
trap ground-state uncertainties being Ax'oApo = hj2. Parameters 
are A 2 = 1CT 5 , 772/771 = 0.75, v = 0.16, A = 0, and S = 25 
(1.000.000 phase-space shifts sampled on a 200x200 grid). Transi- 
tion dipole moments are assumed to be perpendicular to the chosen 
motional axis. 



validity of the central-limit theorem that would predict a per- 
fectly Gaussian profile of the recoil density. However, the 
individual recoil shifts in phase space are not precisely in- 
dependent random variables. The phases, i.e. directions in 
phase space, of subsequent recoils are correlated by means 
of the waiting-time distribution d28l l. The detailed form of 
this distribution seems to be lost in the averaging process, its 
main characteristics, however, a finite mean waiting time be- 
tween subsequent photon emissions, leaves a correlation be- 
tween subsequent recoil shifts. This effect appears as sub- 
stantial deviation from a Gaussian profile in the distribution 
of quadratures shown in Fig. [3(a). In fact, this recoil den- 
sity perfectly approaches a exp(— \x\) function as can be seen 
from the corresponding logarithmic plot in Fig. [3](b). 

It has been shown that fluorescence scattering of a free atom 
generates, after many photon emissions, a Gaussian momen- 
tum distribution j^, @] ■ The above results show that this sit- 
uation corresponds to the limit v — > 0, i.e., when switching 
off the trap potential. Then, the spectral waiting-time distri- 
bution enters Eqs. d73l l and d74T i at zero frequency and equals 
unity due to the normalization of the waiting-time distribution, 
w(0) = 1. Thus, in this particular case the laser parameters 
do not appear, and neither does the waiting-time distribution 
- except its normalization property. Each photon recoil shifts 
the atomic momentum in the momentum's original direction, 
and the length of each shift is a statistically independent ran- 
dom number, whose statistics is given by the dipole radiation 
characteristics. The application of the central-limit theorem 
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Figure 3: (a) Phase-independent distributions of quadratures of p(a) 
for (resonant) fluorescence scattering; parameters are the same as 
in Fig. [2] (1.000.000 phase-space shifts sampled in 200 bins in the 
interval [-1000: 1000]); (b) Logarithmic plot of the same function. 



is here perfectly justified given a sufficiently large number of 
spontaneous photon emissions. 



B. Optical pumping with few photon emissions 

On the other hand, in the opposite limit A2 — > 1, the atom 
spontaneously emits one photon only and reaches its final 
pumped state. In this limit, the joint photon emission prob- 
ability density d32b becomes 



(108) 



lim w n (t n , ...,tx) = S„ A w{ti - t Q ), 



from which the recoil density ( |37T i is derived for completed 
pumping (t — ► 00, to = 0) and making use of d35l ) as 



dt / dsfi2(s)w(t)S[a — ir]2sexp(ii>t)]. (109) 



p(a) 



For high enough laser saturation (S > 1) the waiting-time 
distribution decays at the rate on the order of 7. In an exper- 
iment, the trap frequency is typically much lower, v -C 7. 
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Figure 4: Recoil density p(a) versus position and momentum scaled 
as x/(rji Axo) and p/ (771 Apo), respectively, with trap ground-state 
uncertainties being AzoApo = h/2. Parameters are Ai = 10 _J , 
772/771 = 0.75, v = 0.16, A = 0, and S = 25 (1.000.000 phase- 
space shifts sampled on a 200x200 grid). Transition dipole moments 
are assumed to be perpendicular to the selected motional axis. 



Thus, to very good approximation we may discard in Eq. 
( 11091 ) the comparably slow time dependence of the direction 
of the recoil shift, e lvt w 1, perform the integral over the two- 
dimensional delta function and obtain 



(HO) 



where a = q + ip. Thus the recoil density extends only in di- 
rection of momentum in phase space, reproducing the profile 
of the dipole radiation characteristics. This density distribu- 
tion agrees with the profile associated with a free atom after 
single-photon scattering H]. 

In Fig. |4] we approach this limiting case for the branch- 
ing ratio Ai = 1CP 5 . The almost unidirectional scattering in 
phase space is well observed. Moreover, for the correspond- 
ing quadrature, the distribution of maximum rms fluctuation 
shows a dipole radiation characteristics, see solid curve in Fig. 
[5] The secondary and tertiary lobes in the recoil density of 
Fig. |4]mark delayed spontaneous emission, modulated by the 
laser-driven damped Rabi cycles on the |1) «-> |3) transition. 



C. Laser-saturation dependence of the anisotropy 

The anisotropy can be written as a function of three vari- 
ables, 



A(Ai,w, 0™) = 



(l-Xi)w 



(1 - Aiw) 2 + 4AiWsin 2 (0 m /2) 



(HI) 



13 



120000 



80000 



40000 




quadrature 

Figure 5: Distributions of quadratures of p(a) versus the scaled 
quadrature corresponding to maximum (solid) and minimum 
(dashed) rms quadrature fluctuations. Parameters are those of Fig. 
H (1.000.000 phase -space shifts sampled in 200 bins on the interval 

[-2,2]). 
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Figure 6: Anisotropy A as a function of the laser saturation S for 
resonant pumping (A = 0); parameters are v = 0.16, Ai = 0.75. 



Two of these variables, w and cj) w , depend on the laser param- 
eters, i.e., on detuning and saturation of transition |1) *-» |3). 
For resonant pumping (A = 0), these functions of laser satu- 
ration S are obtained as [see Eqs ( ID 131 ) and ( ID 14M 

w(S) = = =, (112) 



(1 



AV 2 



4>w(S) 



arctan 



i>(S 



S-3i> 2 



(113) 



with v = 2^/7. 

Thus, the anisotropy A(S) depends on laser saturation with 
the branching ratio A2 and the scaled trap frequency v being 
parameters. With Eqs. d 1 12b and d 1 1 3b inserted into flllll ). 
this function is 

A 2 5 



A(S) = 



^J(\ 2 S-3i> 2 f + £> 2 (S + 2-z) 2 f 
Its saturated value is 

lim A(S) = 



A 2 



S — »-oo 



(114) 



(115) 



which vanishes for resonant fluorescence scattering, A2 — > 
0. Likewise, the phase may be written as the saturation- 
dependent function 



tan 4>a{S) 



sin 4> W (S) 



(116) 



cos 4> w (S) — Xiw(S) 

The anisotropy, an example of which is shown in Fig. [6] 
may have a local maximum in form of a typical under-damped 
peak. However, it can be shown that this maximum appears at 
the saturation value 

i> 4 + 5i> 2 + 4 



Srr 



3A 2 



(117) 



only if the condition 



A 2 > 



2-v 2 



(118) 



holds. Otherwise the global maximum is attained at full 
saturation 5 — * 00. For example, given the typical value 
v = 2^/7 = 0.16, a branching ratio A 2 > 0.66 is required to 
observe the maximum anisotropy at extremely high laser sat- 
uration. Typically, this regime is inaccessible in experiment, 
so that usually the anisotropy is a monotonically increasing 
function of the saturation, as shown in Fig. [6] 



D. Optimized optical pumping 

In many applications one intends to minimize the distur- 
bance of optical pumping on an initially given vibrational 
quantum state, in order to only affect the electronic system 
in the desired way. Optimization depends, however, on the 
specific type of disturbance that is supposed to be minimized. 
As we have seen in Sec. IIV Bl the various vibrational mo- 
ments depend quite differently on the laser parameters, some 
even being independent of the kind of laser excitation. Thus, 
we may present optimum values for optical pumping in a com- 
mon case. To discuss the worst-case scenario, we start from an 
atom with initial population in the unpumped state Pi (0) = 1, 

For keeping the rise of the mean vibrational excitation (i.e. 
heating) as small as possible during optical pumping, as in- 
dicated by Eq. j93l , the moment n p , i.e., the mean recoil- 
induced excitation, has to be minimized. This moment does 
not depend on laser parameters, cf. Eq. d68l ), but is a fixed 
value determined by the Lamb-Dicke parameters, the branch- 
ing ratio and the second moment of the dipole radiation char- 
acteristics. It is thus irrelevant whether optical pumping is 
performed resonantly or off-resonantly, or whether a low or 
high saturation is employed. 
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The width of the resulting distribution over vibrational 
quantum numbers, however, does depend on the laser param- 
eters. While the atom is initially in its vibrational ground state 
| n = 0), the final spread of vibrational quantum numbers in 
the pumped state is given by Eq. (1102b as 

([AnM] 2 ) 2 =^ + ^((cos0) 2 )? (119) 

+ (r£C 2 + ^fCi) +2n p ^77i((cos#) 2 ) 1 Acos0 j4 

where we have defined 

Ci, 2 = ((cos#) 4 ) lj2 - ((cos9) 2 )l 2 . (120) 

This is minimized for A cos 4>a being as negative as possible. 
For resonant pumping (A = 0) the anisotropy parameters are 
only functions of the saturation S. Thus, we may search for 
the minimum of A(S) cos <Pa(S). Extrema of this function 
follow from the condition 

A'(S) 

-jrd = <I>'a(S) tan MS), (121) 

where the primes indicate derivatives. The equation can be 
shown to have only negative solutions, whereas the saturation 
by definition is a semi-positive number. Since the global min- 
imum is located at S = 0, where A cos 4>a = 0, a minimum 
spread of the vibrational numbers is obtained when saturation 
is kept as low as possible. Thus, slow optical pumping is re- 
quired for minimum width of the vibrational number distribu- 
tion. Optimization with respect to other moments, however, 
requires other values of saturation. 

VI. MOTIVATION AND CONTEXT 
A. The center-of-mass effects in optical pumping 

Optical pumping is a time-honoured subject. A wealth of 
experiments have required various approaches of modelling. 
However, these theories usually limited their models to the 
description of the atom's electronic dynamics driven by the 
light. The center-of-mass dynamics, modified by the unavoid- 
able radiative recoil, was usually neglected on the basis of the 
assumption that atomic collisions and in particular the spon- 
taneous emission intrinsically tied to the pumping process 
would average an atom's state to an unspectacular slightly 
broadened Gaussian distribution in phase space that might at 
most render energy levels and transitions broadened by an al- 
most unobservable amount. 

The results of the above calculations show that this atti- 
tude is an oversimplification. For complete description of OP, 
the demonstrated modelling of the recoil effects is indispens- 
able. These effects leave, in general, a pumped atom in a non- 
Gaussian state with an anisotropic distribution in phase space. 
The model reveals a continuous transit from maximum asym- 
metry to full symmetry going along with the number of reiter- 
ated scattering events increasing up to the final pumping event. 



Consequently, the center-of-mass effect of optical pumping is 
linked in a natural way to that of resonance scattering on the 
pump-excited line: The latter phenomenon is described by an 
infinite series of scattering events, lacking the final sponta- 
neous decay into the pumped state. When the pump light is 
detuned off resonance, the accumulated recoil turns out to be 
a detrimental heating contribution to the effect of laser side- 
band cooling of the atom weakly bound in a trap. (The wanted 
contribution is from the recoil by the pump light.) For an atom 
strongly bound to the trap, however, the momentum of recoil 
is tranferred to the entire trap. This process is analogous with 
the origin of narrow-band fluorescence emission in the Mtiss- 
bauer effect. 



B. Optical pumping of vapors 

Optical pumping of trapped or free atomic vapors can be de- 
scribed by our model, except for ultra low temperature when 
the gas is degenerate. Elastic collisions will then generate 
scattering among vibrational levels, to effect thermalisation 
of the vapor. 

The atoms will be in a non-equilibrium state after opti- 
cal pumping. As with the electronic inversion, this also ap- 
plies to the vibrational state. Both the particular shape and 
the anisotropy of the density of recoil shifts is transmitted in 
optical pumping to the final vibrational phase-space distribu- 
tion. Such a distribution, in general, is far from a thermal 
distribution corresponding to equilibrium. In particular, the 
anisotropy in phase space indicates that kinetic and potential 
energies of the ensemble are unbalanced, which manifests it- 
self by the atomic cloud breathing at twice the trap frequency. 

Atomic collisions redistribute populations in the vibrational 
levels and thus rebalance kinetic versus potential energy to ar- 
rive at a thermalized equilibrium state. This process happens 
on a time scale that is determined by the atomic scattering 
cross section. Given that this time scale is large enough to be 
detected experimentally, the observation of the cloud should 
reveal damping of the cloud's breathing. 



C. Observability of anisotropy and non-Gaussian shape in 
phase space 

The prominent features of the results presented here are the 
anisotropy and the non-Gaussian shape of the density of recoil 
shifts in phase space. They can be directly observed by using 
a simplified version of a method, proposed for the reconstruc- 
tion of the vibrational quantum state l46l |47||. that has been 
implemented experimentally in recent years 114811 . 

After optical pumping a ground-state atom whose ini- 
tial phase-space distribution is rotationally symmetric (e.g. 
any vibrational number state or incoherent mixtures of these 
states), the |1) «-> |2) transition is bichromatically driven in 
stimulated Raman configuration on the first red and the first 
blue vibrational sideband. Application of this laser light for a 
finite and adjustable time duration t maps information on the 
distribution of quadratures onto the occupation of the ground 
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state |1). Given a laser phase shift <fi introduced in one of 
the two Raman transitions, the ground state occupation reads 



Pi{tA) = \- \ I 



(122) 



where the distribution p(q, (j>) is the probability density of ob- 
serving the quadrature q at phase <fi, as shown in Figs [3] and [5] 
in Sec. H 

The occupation in the ground state Pi is obtained by prob- 
ing the atom for fluorescence of the transition from |1) to an 
auxiliary level. Thus, varying the time duration r and Fourier 
transforming the measured time-dependent ground state occu- 
pation, all distributions of the phase-dependent quadrature are 
obtained. The appearance of phase-dependent distributions of 
quadratures would prove the predicted anisotropy of the den- 
sity of recoil shifts. In addition, the shape of the measured 
distributions could be tested for the deviation from a Gaussian 
profile. 



VII. SUMMARY AND CONCLUSIONS 

In this paper it was shown how recoil effects that naturally 
appear with optical pumping are analytically modelled for a 
harmonically bound atom. Due to the trail of spontaneously 
emitted photons being finite, a quantum trajectory-type of ap- 
proach is well adapted to this analytic procedure. All effects 
on the quantum motion of the atom depend on a single func- 
tion, the recoil density in phase space. Its moments are ob- 
tained analytically from the dipole radiation characteristics 
and the waiting-time distribution of the atom's light scatter- 
ing. 

It was shown that even approaching the case of fluorescence 
scattering the recoil density does not become a Gaussian func- 
tion in the motional phase space, as one would expect when 
taking resort to the central-limit theorem. Instead, it is a distri- 
bution of the form exp(— \x\), that emerges from correlations 
of emission times and the concomitant correlations of the re- 
coil directions in phase space of subsequent spontaneous pho- 
ton emissions. In the opposite extreme, when only few pho- 
tons are scattered during the full pump cycle, the recoil density 
is characterized by the spatial distribution of dipole radiation. 

In general, the recoil density is not isotropically distributed 
but reveals a pronounced directional structure. For resonant 
pumping, this anisotropy can be maximized adjusting the laser 
saturation - in most relevant cases by increasing it. 

Optimum laser saturation and detuning for minimum detri- 
mental heating depend on the particular property that is to be 
protected from the recoil effects. Whereas the mean increase 
of vibrational excitation cannot be varied, the rms fluctuations 
of the vibrational number can be affected by the laser param- 
eters, as is the case also for the rms fluctuations of the phase- 
dependent quadrature. However, optimum laser parameters 
depend on the initial motional quantum state and on the par- 
ticular phase of a quadrature to become protected. There ex- 
ists a complex interplay among the initial motional quantum 



state, the particular fluctuations to be minimized, and the laser 
parameters that fit in with the given requirements. 
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Appendix A: DERIVATIVES OF THE CHARACTERISTIC 
FUNCTION FOR INSTANTANEOUS RECOIL SHIFTS 

The characteristic function associated with the probability 
density of instantaneous recoil shifts is 



P a (a,t) = fj^ [2?7 a |a|cos(j/t- arga)] 



where v is the vibrational frequency. Its (fc, ^)-fold derivatives 
are 



d V / d x 



a=0 



(Vae^) 1 ( Va e- ivt ) k t£ +l H0), 



a=0 



where 



^+0(0) = H) 



k+l 



dsfi a (s)s 



k+l 



H) fc+i ((cos#) fe +') a . 



so that 



a=0 



x ((cos#) fc+i ) a .(Al) 



The dipole radiation characteristics is obtained from Eq. © 



as 



cos 2 a 



} 



(A2) 



where a is the angle between n a and n(£l) and thus 

cos a = cos^cos^a + sin#sin# a cos(0 — <fi a ), (A3) 

where 6 a and (j> a are the spherical angles of the transient 
dipole moment = (3|d|a). Inserting Eq. ( IA3I ) into ( IA2b 
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and performing the integration the dipole characteristics re- 
sults as 



The rms spread of photon numbers is obtained via the 
second-order moment 



M») = r [1 + cos 2 9 a + s 2 (1 - 3 cos 2 e a )} . (A4) 
o 

It is thus a symmetric function, fj, a (—s) = /U (s), that depends 
on the angle 9 a between the transition dipole moment and the 
direction of motion. 



Appendix B: PHOTON-COUNTING STATISTICS 

From the joint probability of spontaneous photon emis- 
sions, the photon-counting statistics derives from the proba- 
bility for n spontaneous emissions during the time interval 

Mo], 

pt pt2 

P„(t,to)= / dt n ... dtiw n (t n , . . . (n>0) 

J to J to 

(Bl) 

and Po(t, to) = 1 — X]^=i Pn(t, to). The photon-counting 
statistics for a complete optical-pumping process, P n , is ob- 
tained in the limit P n = lim^oo P(t, to) which reads 

POO 

P n = dt n ... dt x w n (t n , ...,ti) 

J to J to 

dtiW (t n - tn-i) ■ ■ ■ 



A 2 A™ 1 / dt r , 



to J to 

w(t 1 - 1 ), 



(B2) 



where Eq. (132b has been used. The waiting-time distributions 
are defined for positive time arguments only, so that we set 
w(t) = for t < and extend the upper integration limits 
to infinity. With the substitution of integration variables = 
tk — tfc-i (k = 1, . . . , n) the photon-count statistics writes 



Pn = AaA?- 1 / dr n 



= A 2 A™ _1 (n>0), 



driw(T n ) . ..to(n) 



and Po = 0. Here, the unit normalization of the waiting-time 
distribution has been used. 

The average number of spontaneously emitted photons is 
therefore 



00 



(" P h> = J2 nPn = J2 nX2X i~ 1 > (B3) 

n— n— 1 



which can be expressed as the derivative 

n— 1 

so that the average photon number is 

, 1 



A 2 



(B4) 



(n 2 



pl 



n=l 

71— \ 

n=l 

A2 aA7 ~ V T^aT 

d f d \\ Ai 



J 9Ai V^Ai 
9/9 



= A 2 
= A 2 



9Ai V9Ai 1-Ai 1-Ai 
d Ai 



9Ai (1-Ai) 2 
1 + Ai 



Ai 



Therefore, the rms spread reads 



(An! 



phi \2' 



Its relative value reads 



( A «ph)/<"ph) = T^" 
A 2 

which indicates sub- or super-Poissonian photon statistics, de- 
pendent on the branching ratios. 



Appendix C: MOMENTS OF THE DENSITY OF RECOIL 
SHIFTS 



The function to be evaluated is given by 



(a* k a l ) p = (-l) l J2 / *„.../ dt lWn {t n , . . . M) 

n—1 

< {d*dl* p 2 {a,t n )p 1 (a 1 t n _ 1 )...p 1 {a,t n - 1 ) } . (CI) 



1. Moment (a 2l ) p 



From identity dCll ), the moment (a 2l ) p is given by the ex- 
pression 



00 



)p = X) / dt n ... dtiw n {t n , . . . ,t{) (C2) 

n=l ^ *° •'to 

x |9 2 i [p 2 (a ! i n )p 1 (Q;,t„_i)...p 1 (a ! i„_i) } a _ Q - 
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The non-vanishing second derivative has to be applied to all n 
factors resulting in n terms again. The first term is 



and the normalization of the waiting-time distribution, the n- 
fold integral therefore results as 



E / dt n ... dtiui„(t„, . . . ,ti) 

n=l"'*o J to 



x | [|d 2 ip 2 (a,t„) ^(ajin-i) . . .^(a,^-!) | _ , 

where the derivative is (— irj2) 21 {{cos9) 2l )2& 2llvtrl , cf. Eq. 
(I Alt , so that it reads 



xw n (t nj ...,ti)e au,rt ». (C3) 



The sum of the other n — 1 terms is 



E / dt n ... dtiw n (t n , . . . , tx) |p (oi, t ra ) 

x p (a,i„_i) . . . 9^,p (a,t p ) . . .p (a, i n -i) ! 

— L — - 1 - J — - 1 - J a— 

oo n— 1 />oc />i2 

= (-*7 ?1 ) 2i ((cos0) 2 ') 1 EE/ dt n ... dh 

n=lp=l"'*o •'*0 



x«7„(t„,...,t 1 )e 2ii ^, (C4) 

where the derivative was identified as 
{-irixf 1 {{cos e) 2l )ie 2Uvt K 
In both expressions (IC3t and (IC4b the ra-fold integral 



f/f. 



dt 1 w n (t n ,...,t p ,... 1 t 1 )e ZlWtp 



t 



with p = 1, . . . , n appears. Inserting the definition of the joint 
probability density for n spontaneous emissions d32l i. we find 



POO 



dt„ 



.. / dtiw n (t n ,...,t p ,...,t 1 )e* a '** 

Jto 



J to JtQ 

= A 2 A r 1 l - 1 e 2lZl/t0 \w{2lv)] p . (C7) 
When Eq. dC7| > is inserted into the first term dC3l ), this term 



writes 



\\W_ 



, oo 

Hr, 2 ) 2 '((cos^) 2 ') 2 e 2 ^^V[A : 

•*i — ; 

n—l 

= H» 2 ) 2 ' ((cos e) 2 ') 2 e 2 ^° 1 ^ (2M 



l{2lv)\ n 



The sum of the other 



n — l terms 



\ x w{2lv)' 
[Eq. ( IC4H becomes 



H«if((< 



OO 

cosei) 2 ')^ 2 ^* 



n-1 



oo 

Ew _1 E 

n=l 



= (-*7i) a, <(< 



E H2M] p 

0=1 

, OO 

cos9) 2l ) 1 e 2a ^^ V 



1,-1 p— 1 

A 2 ^ - |ta(2M]' 

Ai ^ 1 1 - w{2lv) 

n—l — v ' 

(cos^) ie ^ /f 2 ;;;. . 

1 — \xw_\2lv) 
Thus, the moment ( IC2t of the density of recoil shifts is 



= (-%f((c 



(a 



2 '\ 



(-1)' ^((c OS 0) 2i ) 2 + ^f((cos^) 2 ')^ 



A 2 u> 
x — 



.J2lv) c2ilut( 
\iw{2lv) 



(C8) 



A 2 A 



2*1 



dt„... / dii™(i n -f n _i)...w(ii-i K l(,/t *\ 

Jto 

(C5) 

In the vein of Eq. ( 1B21 i. the upper limits of the integrals are 
set to infinity. With substitution of the n integration variables, 
r n = t n — i„_i, expression ( !C5t gives 



X2X n,-l e 2iluto / 

Jo 

xe 2ilu(r p + ...+T 1 ) 



driw{T n ) . ..w{n) 



2. Moment (a* 2 a 2 ) p 
From Eq. ( ICU the moment {a* 2 a 2 ) p has the form 



{a* 2 a 2 



00 pOO ft 

p = E / • ■ ' / 



dtiw n {t n , . . . ,ti) 



where t p = r„ + r„_i 



ri + to is used. Given the 



definition of the spectral waiting-time distribution 



w{lo) = / dTw{r)e luJT , 
Jo 



(C6) 



The non-vanishing derivatives acting on a single factor 
P a (a, t) are those of second and fourth order, i.e. the deriva- 
tives 9 2 .p , <9 2 p , d a d a *p , and <9 2 i9 2 »p . Therefore, the 
moment can be written as 
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/ *2„2\ 

(a a ) 





rt2 




/ dt n . . 


. / dtiw n {t n , . 


...tl){ 


/to 


Jto 





EE {[ 5a * 5a £a,, da'daP^Jfititq) + dlp^^ (a, t p ) ^.p^ (a, t 9 ) J > , (C9) 



0=1 q^p 



a=0 



where a p = 1 forp = 1, . . . , n — 1 and a p = 2 forp = n. The 
fourth-order derivatives that corresponds to the first term in 
the braces, produce terms such as 7/ 4 ((cos6>) 4 ) a . They occur 
once with a = 2 and n — 1 times with a = 1, to yield, after 
time integrations, the contribution 

% 4 ((cos0) 4 ) 2 + ^^((cos0) 4 ) 1 . 

The next term in the braces of Eq. ( !C9t results as the time- 
independent double sum 

EE<<« cos0 ) 2 U( cos0 ) 2 >^ 

p=l q^p 

= 2r)Wi(n- I)<(cos0) 2 > 2 ((cos0) 2 )i 
+7 ? 4 (7i-l)(n-2)((cos0) 2 ) 2 , 

which, after time integrations in Eq. flC91 l, adds to the moment 
the contribution 

2^1 «ftph) ~ 1) ((cos0) 2 ) 2 ((cos0) 2 ) 1 

+ [r, 4 «n p h> - 1) ((V) - 2) + ^ 4 <[An ph ] 2 )] <(cos#) 2 ) 2 

= 2r / 2 ^((cos0) 2 ) 1 (% 2 ((costf) 2 ) 2 + ^((costf) 2 )^ . 

The last term in the curved bracket of Eq. ( !C9t produces 
the double sum 

EE«((co S «) 2 ),((cos«) 2 )„,e & ^y 

p=X q^ p 

n p—1 

= 25R^^ ?7 2 p?7 2 ((cos0) 2 ) Qp ((co S 0) 2 ) 1 e 2 ^^), 

p=l q=\ 

which generates the contribution 

oo n p—1 

m E E E « cos ) 2 M( cos ^) 2 ) i ( C1 °) 

ri— 1 p—1 q—1 

x dtn... dt 1 w n {t ni ... 1 t 1 )e 2lv(tp ~ tq) - 

Jto Jto 



The time integrals, 



dt, 



to 



' dtiw n (t n ,...,t 1 )e 3iv <- t >- t *\ 



with p > q, can be transformed as above into 

/•oo />oo 

AaA?- 1 / dr n ...\ dT 1 w{r n )...w{Tx)e 2iv ^ + - +T "+^ 
Jo Jo 

Thus, the contribution dClOb becomes 

00 n p—1 

m E E E < ? ?i « cos ^) 2 M( cos ^) 2 ) 1 

n— 1 p—1 g— 1 

xA 2 A™- 1 [ m (2 J /)]^ 9 . 

Using the substitution g — > q' = p — q and performing the 
sum over q', we have 

00 n 

^ E E <^ 2 « cos °) 2 )^ « cos ^) 2 ) 1 ^Ar 1 

n— 1 p—1 

w{2v) 



1 - u;(2z/) 



{l-t^r 1 }. 



In the sum over p we separate the term with p = n from the 
terms p = 1, . . . , n — 1 and obtain 

2^£ ?7 2 ((co Se ) 2 ) 1 A 2 Ar- T ^- y 

n=l — ^ ' 

xj^Ucos^jl-^^r 1 } 

1 - [w(2i/)]™- 1 1 1 



+^((cos^)iUn-l) 



The sum over n produces terms like 



1 - w{2v) 



OO . 

E^cAixr-^ 

n=l 



1 — XiX 



where x may be 1 or w_(2v). Thus, after performing this sum, 
we arrive at 
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25Rt7?((cos6») ; 
y.Ll{(cose)% 



w{2u) 
l-w{2u) 



1 - 



1 - \iw{2v) 



The brackets can be simplified 
[l ^ 



to 



= Xi[l-w(2v)] 
1 - \iw(2v)\ 1 - \iw(2v) ' 

and (h p h) = I/A2, see Eq. (IB4t . so that we obtain 
contribution ( IC1 U 



for the 



2^((cos 



1 - Aiw(2z/)_ 
The complete results is therefore 



df)i (% 2 ((cosfl) 2 ) 2 + ^ 2 ((cos0) 2 



(a* 2 « 2 ) P = 



??2 4 ((COS Of) 2 + ((COS 0) 4 >! 

A 2 



+2rA(co S fl) 2 ) 1 
A 2 



^ 2 2 ((co S 0) 2 ) 2 + ^7 ? 2 ((cos0) 2 



1 + 9? 



A 2 w(2;/) 
1 - Aiu;(2i/) 



(C12) 



Appendix D: SPECTRAL WAITING-TIME DISTRIBUTION 
FOR RESONANT PUMPING 



The waiting-time distribution w(i) is defined by 



«>(*) = (7i +7a) (3|t>cff(*)|l) 



(Dl) 



[Eq. (I28H, The time evolution U c g(t) is governed by the non- 
Hermitean Hamiltonian 



H oS = -(«o-i,+ + K*a-i l _)+nA|l)(l| 
ih \ -> 



(D2) 



[Eq. (9)]. The transition amplitude |1) — > |3) is the projection 
of the solution of the equation of motion 



^ 7^— = i/eff|V'W): 



(D3) 



on state |3) with the initial condition \ip(t = 0)) = |1). Eq. 
is decomposed in the set of coupled equations for the 



+ ^((Cosfl) 2 )! ((h ph ) - 1) - -— 

L 1 — w(2i/J 



1 - 



A 2 



1 - Xiw(2u)_ 



(Cll) 



components ipi(t) = (i\ip(t)), 



Mt) = -—i> 3 (t)-iAMt), (D4) 



Mt) = 0, 

Mi) = -yV'iW -7V»3(t). 



(D5) 
(D6) 



where 7 = (71 +7 2 )/2. The second derivative $3(4) therefore 
results as 

which, with the use of Eq. ( ID4I ). writes in the resonant case 
(A = 0) 



d 2 t + 7 a t + 



Mt) = 0. 



(D7) 



With the Ansatz 

ip3(t) = <f> 3 {t)cxp 
the differential equation for cfe (t) is obtained: 



2 



hit) = 

The general solution is 

Mt) 

with 





K 


2' 


[(I) 2 ' 


2 





»(«). 



6e 



-A/ 



(D8) 



(D9) 



(D10) 



where S = \k\ 2 /j 2 is the laser-saturation parameter. 
Thus, the time-evolved amplitude of state |3) reads 



Tp 3 (t) = (ae xt +be- xt )e-^ 2 . 

With the initial conditions ^3(0) = and ipi(Q) = 1, one 
obtains, via Eq. dD61 >. ^3(0) = —in/2, so that one finds the 
constants of integration, a = — b = — i/t/(4A). Including 
these values, the above probability amplitude is 

Mt) - J[_ s sinh (| VT^S^j e^' 2 . (Dll) 
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The waiting-time distribution w(t) = 2j\ip3(t)\ 2 is then given 
by 



J{t) = TT^sf 



it 



sinh ( yVT^5 



-it 



(D12) 



Note that this distribution is normalized to unity: 

/ dtw(t) = 1. 
Jo 

The spectral waiting-time distribution results from the 
Fourier transform of Eq. dD12b . 

/>oo 

w{u) = / dtw{t)e lut . 
Jo 

It is obtained by straightforward calculation as 

7 3 S 



w(uj) = 



( 7 -^)[ 7 2 (S-l) + ( 7 -iw) 2 



1, S<1, 
(-1), S> 1. 



Thus, the required modulus and phase at twice the trap fre- 
quency become 



5 w (sat) 

- ff)1 3) 

V(5 - l) 2 + 2(5 - 1)(1 - z> 2 ) + (1 + z> 2 ) 2 



tan^ = s _ 3 ~ 2 tan^ at > 
where the saturated values are 



(D14) 



w (sat) = lim w = l/VT + ^2, 

s — >oc 



= lim 



= arctani/, 



with i> = 2v/ r y. 



[1] A. Kastler, J. Phys. 11, 255 (1950). 

[2] Symposium "Alfred Kastler", Ann. Phys. Fr. 10, No 6 (1985). 
[3] P. Kapitza and P.A.M. Dirac, Proc. Cam. Phil. Soc. 29, 297 
(1933). 

[4] A. Ashkin, Phys. Rev. Lett. 24, 156 (1970). 
[5] O. Frisch, Z. Physik 86, 42 (1933). 

[6] A. Yu. Pusep, Zh. Eksp. Teor. Fiz. 70, 851 (1976) [Sov. Phys. 

JETP43, 441 (1976)]. 
[7] J.F Lam and PR. Berman, Phys. Rev. A 14, 1683 (1976). 
[8] L. Mandel, J. Opt. (Paris) 10, 51 (1979). 
[9] V.S. Letokhov and V.G. Minogin, Phys. Rep. 73, 1 (1981). 
[10] T. Hansen and A. Schawlow, Opt. Commun. 13, 68 (1975). 
[11] V.S. Letokhov, V.G. Minogin, and B.D. Pavlik, Opt. Commun. 

19, 72 (1976); ibid. Sov. Phys. JETP 45, 698 (1977). 
[12] S.V. Andreyev, VI. Balykin, V.S. Letokhov, and V.G. Minogin, 

JETP Letters 34, 442 (1981); ibid. JETP 55, 828 (1982). 
[13] S. Chu, Rev. Mod. Phys. 70, 685 (1998); C.N. Cohen- 
Tannoudji, ibid. 70, 707 (1998); W.D. Phillips, ibid. 70, 721 
(1998). 

[14] H. Metcalf and P. van der Straten, Laser Cooling and Trapping 

(Springer Verlag, Berlin, New York, 1999). 
[15] W. Neuhauser, M. Hohenstatt, PE. Toschek, and H.G. Dehmelt, 

Phys. Rev. Lett. 41, 233 (1978). 
[16] D.J. Wineland, R.E. Dnillinger, and FL. Walls, Phys. Rev. Lett. 

40, 1639 (1978). 

[17] W. Neuhauser, M. Hohenstatt, PE. Toschek, and H. Dehmelt, 
Phys. Rev. A 22, 1137 (1980). 

[18] A.D. Boozer, A. Boca, R. Miller, T.E. Northrup, and H.J. Kim- 
ble, Phys. Rev. Lett. 97, 083602 (2006). 

[19] S. Nussmann, K. Murr, M. Hijlkema, B. Weber, A. Kuhn, and 
G. Rempe, Nature Physics 1, 122 (2005). 

[20] P. Meystre, Atom Optics (Springer Verlag, New York, 2001). 

[21] V.S. Letokhov, Laser Control of Atoms and Molecules (Oxford 
University Press, New York, 2007). 

[22] E. Joos, H.D. Zeh et al., Decoherence and the Appearance of the 
Classical World in Quantum Theory (Springer Verlag, Berlin, 
2003). 

[23] R. Blatt, J.I. Cirac, and P. Zoller, Phys. Rev. A 52, 518 (1995). 



[24] S. Wallentowitz, W. Vogel, I. Siemers, and PE. Toschek, Phys. 

Rev. A 54, 943 (1996). 
[25] Th. Beth and G. Leuchs, eds. Quantum Information Processing 

(Wiley- VCH, Weinheim, 2005). 
[26] A.M. Steane, Rep. Prog. Phys. 61, 117 (1998). 
[27] M. Sasura and A.M. Steane, Phys. Rev. A 67, 062318 (2003). 
[28] D. Haubrich, H. Schadwinkel, F. Strauch, B. Ueberholz, R. 

Wynands, and D. Meschede, Europhys. Lett. 34, 663 (1996). 
[29] R. Bonifacio, L. De Salvo, L.M. Narducci, and E.J. D'Angelo, 

Phys. Rev. A 50, 1716 (1994). 
[30] D. Kruse, C. von Cube, C. Zimmermann, and Ph.W. Courteille, 

Phys. Rev. Lett. 91, 183601 (2003). 
[31] GC. Hegerfeldt and T.S. Wilser, in Proceedings of the II. In- 
ternational Wigner Symposium, 1991, eds H.D. Doebner, W. 

Sherer, and F. Schroeck (World Scientific, Singapore, 1992), p. 

104. 

[32] C.W. Gardiner, A.S. Parkins, and P. Zoller, Phys. Rev. A 46, 
4363 (1992). 

[33] J. Dalibard, Y. Castin, and K. M0lmer, Phys. Rev. Lett. 68, 580 
(1992). 

[34] H.J. Carmichael, An Open Systems Approach to Quantum Op- 
tics, Lecture Notes in Physics Vol. ml8 (Springer Verlag, 
Berlin, 1993). 

[35] K. M0lmer, Y. Castin, and J. Dalibard, J. Opt. Soc. Am. B 10, 
524 (1993). 

[36] B.M. Garraway and PL. Knight, Phys. Rev. A 50, 2548 (1994). 
[37] M.B. Plenio and PL. Knight, Rev. Mod. Phys. 70, 101 (1998). 
[38] W. Vogel and S. Wallentowitz, in Coherence and Statistics of 

Photons and Atoms (Wiley, New York, 2001), ed. J. Perina. 
[39] C. Cohen Tannoudji and J. Dalibard, Europhys. Lett. 1, 441 

(1986) . 

[40] P. Zoller, M. Marte, and D.F. Walls, Phys. Rev. A 35, 198 

(1987) . 

[41] M. Porrati and S. Putterman, Phys. Rev. A 39, 3010 (1989). 
[42] M.S. Kim and PL. Knight, Phys. Rev. A 40, 215 (1989). 
[43] K.E. Cahill and R.J. Glauber, Phys. Rev. 177, 1857 (1969). 
[44] K.E. Cahill and R.J. Glauber, Phys. Rev. 177, 1882 (1969). 
[45] C.W. Gardiner, Handbook of stochastic methods (Springer, 



21 



Berlin, 1985), 2nd ed. 
[46] S. Wallentowitz and W. Vogel, Phys. Rev. Lett. 75, 2932 (1995). 
[47] S. Wallentowitz and W. Vogel, Phys. Rev. A 54, 3322 (1996). 
[48] PC. Haljan, K.-A. Brickman, L. Deslauriers, PJ. Lee, and C. 

Monroe, Phys. Rev. Lett. 94, 153602 (2005). 
[49] Without loss of generality, the zero-energy level has been 



shifted here by — h(oJi + 0J2) for simplicity. 
[50] This is a standard expansion, known in quantum optics by the 

term "unraveling" of the master equation. 
[51] This case is impossible to realize, at least for resonant pumping 

(A = 0), as can be seen from the results of App.lDl 



